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We present a new model for a non-topological soliton (NTS) that contains fermions, scalar particles 
and a gauge field. Using a variational approach, we estimate the energy of the localized configuration, 
showing that it can be the lowest energy state of the system for a wide range of parameters. 



(N 
O 
O 
(N 

Oh' 
(D 

m 



> 

in 

cn 
o 



Oh! 
I 

Oh 



X 



PACS numbers: 11.10.Lm,11.15.Tk 



I. INTRODUCTION 



The study of solitons and non-topological solitons 
(NTSs) has a rich history. They have been proposed 
as building blocks for stars and black holes ^, and 
as dark matter candidates |3[ |[ ||. The first models for 
an NTS in 3+1 dimensions were found by Lee and Wick, 
and by Friedberg, Lee and SirHn [|, §. The latter NTS 
contained one real scalar field to set up a false vacuum in 
which a second, complex scalar field was massless, allow- 
ing the NTS to be the lowest energy state for the system 
under certain conditions. 

Coleman and collaborators extended this work to sim- 
pler objects dubbed Q-balls, which contained a single 
complex scalar field that possessed a conserved global 
symmetry ^. More recent work has extended NTSs 
to contain fermions and gauge fields Finally, 
work has been done by Kusenko extending NTSs to some 
supersymmetric field theories, where the corresponding 
false vacuum is set up in the superpotential Q, |5| . 

In this paper we present a new model for an NTS con- 
taining mutually interacting fermions, scalar particles, 
and a U{1) gauge field. The fermions have a Yukawa cou- 
pling to the scalars, and both carry a conserved global 
charge. This model is thus closer to more realistic theo- 
ries of particle physics that fit in or are inspired by ex- 
tensions of the Standard Model. We note that recently a 
fermionic Q-ball model was proposed by Anagnostopou- 
los et al. [|l2j. However, the approach adopted in the 
present manuscript is quite different, as we explicitly 
solve the Dirac equation, as opposed to modeling the 
fermions as a relativistic gas from the outset, and include 
the Yukawa coupling to the scalar field. Our results, nev- 
ertheless, agree with the general conclusions of Ref. [ p^ , 
as we show that, indeed, it is possible to obtain NTSs 
which are the lowest energy state in the system for a wide 
range of parameters. In addition, there is an opportunity 
to extend our work to bring it in line with the nuclear 
bag model, as the expression derived below for the en- 
ergy approximates some models used in nuclear physics. 
Throughout we work in natural units where h = c = 1. 



II. THE NTS 



Consider the Lagrangian 



C = {D,MD^<I^T - -^F^'^F,,, - Um + (1) 



iTp^^ D fitjj — ml 



where (j) is the complex scalar field, = — ieA^ 
is the U(\) covariant derivative, F^i, is the field tensor, 
and J7(|0|) is the potential for the scalar field, ip is the 
4-component spinor, 7p is the Dirac adjoint spinor, 
are the four covariant Dirac matrices, to is a positive 
constant, and is a constant chosen such that when 
\(t>\ = F_ the fermions are massless within the NTS. We 
can express the complex scalar field as two real fields 
using = exp(i6'), to get 



C = \d,fd^f + \f\d,e eA,f U{f) \f,^F'^'^{2) 
-1-1-07^(9^ - ieA^)%l; - to^I - -y^^il). 

We assume the ground state will be spherically symmet- 
ric, and will have no magnetic field and hence no elec- 
tric currents. Therefore, we may choose a gauge where 
Af_i — A(){r). The boundary condition is that ^o(^) 
as r — > cx). In addition, we make the assumption that the 
scalar field oscillates in time with a regular frequency and 
hence 6 = ojt, where w is a positive constant ||ll[. The 
Lagrangian then becomes 



L — An / r dr 



1 



1 



1 



-T:f + Tr^r + 7;f9-U{f) (3) 
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-\-tp^ (w — g)'ilj — ml — 



where g = oj — eAQ(r). The Eulcr-Lagrange equations for 
g, f and "0 a-i'G 



0, 



(4) 
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/" + If + fg' - ^ + = 0, (5) 



h'^id^ - ieA^)ip -mil 



t/j = 0, 



(6) 



where it is understood that the only non-vanishing com- 
ponent of Af^ is Aq. The conserved currents and charges 
are given by 



(7) 



Q = j Jo"°'"t^'a; - 47r y r'^drpg, (8) 



ij f errriion 



(9) 



Therefore, we must choose one of them to equal zero, and 
conventions of a right-handed coordinate system dictate 
that it is the left handed component that must be zero. 
We can now expand in terms of an angular part and a 
radial part by writing [Q 



(j)R = ih(r)VLjir. 



(14) 



where ^jim is a spherical spinor and the indices j,l,m 
are the quantum numbers of total angular momentum, 
orbital angular momentum, and the z-component of 
angular momentum, respectively. In the spherically- 
symmetric ground state we have I = and so j = 
I + 1/2 = 1/2, and the spherical spinor for this case sim- 
plifies to |l| 




(15) 



and we see that indeed our wave-function is spherically 
symmetric. The equations for h{r) become |13[ 



iV = / J'/"™""d3x = An / r^drij^^. (10) 



dh 
dr 



(1 



0, 



(16) 



III. TAKING CARE OF THE FERMIONS 

To proceed, we could take one of two approaches. We 
could use the Fermi gas approach as in ||l^, |l^. Or we 
could attempt to solve the equations of motion directly. 
We attempt the second method here. We begin with the 
Dirac equation for the fermion field. We can write in 
terms of two, 2-component, spin-1/2 spinors in the chiral 
representation as 



(11) 



where 4>r(l) is a right (left)-handed spinor. It is easier to 
proceed if we switch to the non-covariant representation 
of the Dirac equation. Multiplying by = (7°)~^ on 
the left and using that 7* = /3q;% we get the equation 



where V/(r) — eAQ{r) is the potential for the fermions, 
and M = (1 — f / FJ)m is the mass of the fermion. Fol- 
lowing Lee et al. we will assume that /(r) = constant = 
= F inside the NTS ||ll|. We will also assume that 
■0 is of the simple form if{r,t) = ip{r) e'>ip{—iEt) where 
E is the energy of a single fermion. We therefore get the 
new equation 



Eip = {-ia ■ V -I- Vf{r))ip. 



We desire ip to be spherically symmetric in the ground 
state. Since (pR and (pL must have opposite parity un- 
der spatial reflection, they cannot both be symmetric. 



[E-Vj 



0, 



(17) 



where k= — (/ + 1) so here k = — 1 |T^. The solution to 
Eq. ([1^ ) is then simply h — constant = B. However, we 
notice a problem with Eq. (p7|); namely, if V/(r) is non- 
constant then it demands that h — since E must be 
a constant, and the equation must be sastified for all r. 
How are we to get around this problem? We can approx- 
imate the fermionic energy inside the NTS as constant 
by taking the expectation value of the fermionic poten- 
tial inside it. Then we can write E « (E) = (V/(r)) and 
Eq. ( p7| ) may be approximately satisfied inside the NTS. 
Before doing this we note that if h is constant, the expres- 
sion for the conserved charge Eq. (|lO|) becomes simply 



N ^ 



(18) 



(12) 

where R is the radius of the ball. We then get that 



(Vf) = An / 'ip^(73V{r)'iJjr'^dr ^ — j / r'^[uj - g{r)]dr, 
Jo R Jo 

(19) 

where 173 is the third Pauli spin matrix, and our approx- 
imate solution is 



(13) ^{r,t) 



1 /3N 



1/2 



where 6{R — r) is the step function 



exp(^-it{Vf)^0{R-r) 





Vo 



(20) 
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We can now use Eq. (20) in Eq. The solution for 
g{r) is 

( [uj- e'^Q/ATTR]Rsmh{eFr)/rsmh{eFR)+ 
g{r) = <^ r<R 

(21) 

where R is the radius of the soUton. We can now substi- 
tute for pg in Eq. (||) to get 



(22) 



Using our solution for g{r), Eq. ([2l|), we get 



AttR 



X — tanh(x) 



e^A^ tanh(a;) 
4:TrR{x — tanh(a;)) ' 



where x = eFR. 



IV. ESTIMATING THE GAUGED FERMIONIC 
Q-BALL ENERGY 

The gauge invariant energy can be written as 



E = in I r^dr 



2^ +2^^ + 2^ + 



where Ep is the relativistic Fermi energy given by 
Ep = (37r/4)(3/27r)2/3(iV4/3/A). Ignoring surface terms 
0{R^), performing some partial integrations, and using 
our solution for we get, for r < R, 



E < -ujC 
- 2 



-ttR^U{F) + Ep 



We can use Eq. ( p3[ ) to get 



E < 



SttR 



tanh x 



R 



SttR 



-ttR^U{F) - 

tanhx 
X — tanh X 



(25) 



(26) 



where Ci = 37r/4(3/27r)2/3. The next step is to minimize 
this expression with respect to the various parameters 
at fixed TV and Q. We examine the case of small e and 
expand to order e^, where the Laurent series are given by 
x/{x — tanhx) ~ 3/a;^ + 6/5 and tanha;/(x — tanhx) ~ 
3/a;^ + 1/5. We then get an approximate form for E to 
0{e^) of 



E < 



3Q2 



SttF^R^ 



(27) 



3e2g2 3QN e^QN 
20nR ~ 8nF^R3 ^ iOnR ' 



These terms now have an easy physical interpretation. 
The first term is the zero-point energy of the scalar par- 
ticles. The second term is the vacuum volume energy 
of the bag. The third term is the Fermi energy. The 
fourth term is the Coulomb repulsion of the scalar par- 
ticles. The fifth and six terms represent the interactions 
of the fermions with the scalars and gauge fields, which 
may significantly alter the NTS's energy. 

If we minimize this expression with respect to R we 
obtain, writing R^ 



AtiU{F)v^ - V 



= V, 

3e^ 
207r 



87rF2 



{Q-N)^ 0. 



(23) 

The formal solution to this equation is given by 



A = 



1 fiB /3\ 27rn 

- arccos — -\/ -r H 

3 \2A \ A) 3 

n = 0,±1,±2. 

z = y/a, a 



4:nU{F) IGOn^UiF) 
9Q 



B = 



32tt^F'^U(F) 



7(4173), 
(6Q-iV), 

(Q-N). 



(29) 

(30) 
(31) 



Unfortunately, since we must minimize the energy with 
respect to both R and F the details of selecting the cor- 
rect root (or even if a positive real root exists) depend 
on the nature of U{F). Hence, it is not possible to write 
down a general procedure for selecting a root. Once a 
suitable potential has been chosen however, the solution 
follows in a straightforward manner. To illustrate this, 
we solve for the potential 



U{f) 



11 

4 



(32) 



We choose A = 0.444, ^ = .25, e = .1, Q = 10000, 
and M^i, — 2 (we must choose > 3/16 to insure 
that [/(/) > for all / 0), where fi and have 
dimensions of mass {M^ is the mass of the fermions 
in the true vacuum). The energy scale is set by F, 
which has dimensions of mass. We first solve the prob- 
lem with no fermions present, which is identical to the 
scenario considered in p4| |. This gives us an idea of 
a possible range of values for R and F (based on this 
analysis, we set F = 0.48 throughout) to search for 
fermionic NTSs. Writing the energy of the free scalars 
and fermions as -Efrco = mQ + M^N ^ the NTS is stable 
whenever Ej^t s / E^cc < 1. In Fig. |^, we show the ratio 
Ents / Ehcc as a function of increasing fermion number, 
for several values of the radius. From a more detailed 
analysis we can show that, for the parameters used in 
Fig. ^ the condition for the existence of NTSs is satis- 
fied for 35 < R < 111. Clearly, the same sort of range 
search can be performed for any set of paramenters. 
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In Fig. 1^ we show the ratio E^ts / E^cc as a function 
of the NTS radius for different values of the fermionic 
charge N. Clearly, there is a wide range of values for N 
wherein the NTS is the preferred energy configuration. 
For large radii (larger than the range of values displayed 
in the Figure), the NTS energy is independent of N, as 
can be easily seen from Eq. (UTh- 

We conclude with a few remarks about our NTS. First, 
we see that if we use a potential that satisfies Coleman's 
condition |^ min[2t/ /|0^|] < /i^ (/i is the mass of the free 
scalar particles) , in order to set up a false vacuum where 
the fermions can be massless, then we can always find pa- 
rameters where Emts < E^ee- Second, we see that the 
presence of both fermions and the gauge field increases 
the energy and the radius of the NTS, while the attrac- 
tive Yukawa coupling between the fermions and scalars 
decreases both in relation to the ungauged scalar Q-balls 
studied by Coleman and collaborators H, ||. Third, we 
note that the asymptotic form of the energy is 



FIG. 1: Ratio of NTS energy to the energy of free fermions 
and scalars as a function of fermionic charge N. The labels 
indicate the different values of the NTS radius. 



lim E 

R^oo 



;TtR^U{F). 



(33) 
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FIG. 2: Ratio of NTS energy to the energy of free fermions 
and scalars as a function of radius R. The labels indicate the 
different values of the fermionic charge A^. 



Hence, if we scale U{F) in such a way that R^U{F) — 
constant < jEh-cc, our NTS can be stable at arbitrarily 
large radii. Physically, we have obtained a state of mat- 
ter with scalar particles and fermions uniform through- 
out. We thus see that it is possible to form a NTS out of 
gauged fermions and scalar particles, which can be the 
preferred energy state of the system for a wide range of 
parameters. It would be interesting to investigate the 
solutions to the set of coupled equations numerically, ob- 
taining a more detailed analysis of the allowed parameter 
space. 
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